Spontaneous alloying in binary metal microclusters 
- A molecular dynamics study - 
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Microcanonical molecular dynamics study of spontaneous alloying(SA), which is a manifestation 
of fast atomic diffusion in a nano-sized metal cluster, is done in terms of a simple two dimensional 
binary Morse model. Important features observed by Yasuda and Mori are well reproduced in our 
simulation. The temperature dependence and size dependence of SA phenomena are extensively 
explored by examining long time dynamics. The dominant role of negative heat of solution in 
completing S A is also discussed. We point out that a presence of melting surface induces the diffusion 
of core atoms even if they are solid-like. In other words, the surface melting at substantially low 
temperature plays a key role in attaining SA. 

PACS numbers: 05.60.-k, 36.40.Sx, 61.46.4-w, 66.30.Jt, 67.80.Mg 
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I. INTRODUCTION 

Microclusters exhibit neither the properties of bulk nor 
those of molecules. One of the most important features 
of atomic and molecular clusters is the presence of the 
large portions of surface atoms, which provide large fluc- 
tuation in their motion. A lot of interesting static and 
dynamical properties of nano-sized clusters are found by 
many authors during the last two decades 0. If we re- 
strict ourselves to dynamical aspects of microclusters, for 
instance, it was shown that small metal clusters fluctu- 
ate between different multiply twinned and single-crystal 
structures rather than having fixed structures 0. From 
the viewpoint of equilibrium thermodynamics, Ajayan 
pointed out that this is a manifestation of a quasi-molten 
state where the Gibbs free energy surface as a function 
of the cluster morphology is quite shallow ||| . From the 
viewpoint of a molecular dynamics (MD) study, on the 
other hand, Sawada and Sugano showed that the struc- 
tural change of Au clusters observed in experiments is 
regarded as a floppy motion between local minima of a 
potential surface due to the dynamical nature of clus- 
ters [Q. In both cases we may say small clusters suffer 
from anomalously large dynamical fluctuations. Owing 
to the presence of such large fluctuations, as pointed out 
by Sugano, it is hard to give a clear answer for the follow- 
ing naive questions; Are microclusters like solids where 



atoms are oscillating around their respective equilibrium 
positions? Are they like liquids where atoms move dif- 
fusively? Or, are they fluctuating between different solid 
phases during the course of their motion? Q In fact, 
according to the works by Berry and his coworkers in 
their microcanonical MD study of a Lennard-Jones clus- 
ter, there exists the intermediate phase called 'co-existing 
phase' of liquid and solid j|] . While numerical searches of 
stable structures of small microclusters have been done 
extensively, research is left untouched on extremely long 
time dynamics beyond micro seconds, which is responsi- 
ble for diffusion process in clusters. In the present pa- 
per, we discuss a novel fast diffusion process which was 
experimentally discovered by Yasuda and Mori(YM) in 
nano-sized binary metal clusters, because it is a mani- 
festation of an anomalous diffusion process peculiar to a 
microcluster || . The aim of the present work is to real- 
ize spontaneous alloying(SA) with a concise model and 
to elucidate what kind of dynamics dominates SA with 
an extensive numerical study. 

The present paper is organized as follows. In the next 
section we mention the experimental results of SA and 
an outline of our motivation. In Sec. Ill our model for 
the MD simulation of SA is introduced with a physical 
assumption which we made to prepare appropriate ini- 
tial configurations. Numerical results and observations of 
the MD simulation are presented in Sec. IV. Some useful 
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quantities are introduced to characterize atomic fluctua- 
tion and rearrangement in a cluster. A comparison be- 
tween experimental and numerical results are discussed in 
Sec.V. In Sec. VI the differences between surface and core 
atoms in a cluster are emphasized by paying attention to 
the activity of atoms manifested by their fluctuation and 
rearrangement. A special emphasis is put on the impor- 
tant role of surface melting in SA. Lastly, in Sec. VII we 
briefly make a concluding statement regarding the inter- 
pretation of the results. 



II. SUMMARY OF YM'S EXPERIMENTAL 
RESULTS AND AN UNUSUAL FEATURE OF SA 

A. YM's experimental results 

In 1992 a novel transport phenomenon in a nano- 
sized metal alloy cluster was reported by YM [^J. By 
using an evaporator they deposited individual solute 
atoms (copper) on the surface of host nano-sized clus- 
ters on amorphous carbon film at room temperature and 
observed the alloying behavior with a transmission elec- 
tron microscope. In Fig. 1(a) their in situ observation is 
schematically described. In j|] they demonstrated that 
gold clusters promptly changed into highly concentrated, 
homogeneously mixed (Au — Cu) alloy clusters. This pro- 
cess is termed as spontaneous alloying(SA). SA is simi- 
larly observed in many nano-sized binary clusters, such 
as (Au-Ni),(In-Sb),(Au-Zn), and (Au-Al) §. They 
examined the presence and absence of SA for clusters of 
various sizes. 
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in bulk crystalline alloys. By use of the simple rela- 
tion, x = V Dt, between the diffusion coefficient D and 
the time t, needed to achieve diffusion of solute atoms 
across the distance x, YM roughly evaluated the value 
1.1 x 10 -19 [m 2 /s] as a lower limit. Note that the diffu- 
sion constant of copper in the bulk gold is known to be 
2.4 x 10- 28 [m 2 /s] at 300[K] §. 

(2) Negative heat of solution plays an important role in 
enhancing and supressing the SA process. Indeed, SA 
has never been observed in the combination of the solute 
and host atoms with sufficiently large positive heat of so- 
lution. However, it is worth mentioning that SA occurs 
even with the species with almost null heat of solution. 

(3) Temperature is also an important factor control- 
ling SA. At relatively high temperature(T ~ 245 [K]), 
Cu atoms can dissolve well into the center of a 4nm- 
sized Au cluster, whereas at medium temperatures(T ~ 
215 [K], 165 [K]), the dissolution of copper takes place only 
over a limited, shell-shaped region beneath the surface 
of a 4-nm sized cluster, and the thickness of that region 
where a solid solution is formed decreases with the de- 
crease in temperature. 

(4) With increase in cluster size the occurrence of rapid 
SA is suppressed. In Au clusters of approximately lOnm 
in the mean size, rapid alloying of Cu takes place only at 
the shell-shaped region beneath the free surface of an in- 
dividual cluster and pure gold was retained at the central 
region of the cluster. In Au clusters of approximately 
30nm in the mean size, no rapid alloying of Cu does not 
take place. It should be stressed that the critical size of 
the SA increases with the negative heat of solution and 
temperature. 

(5) SA takes place in a solid phase, which was confirmed 
by the fact that no changes were observed in the multiply 
twinned structure of the host cluster during the alloying 
process. Although the electron beam heating seemingly 
brings considerable influence on SA, the estimated mag- 
nitude of the temperature rise in atom clusters was of 
the order of 10 [If], which causes no significant effect on 
SA. 



B. What is unusual in SA? 



Fig. 1 A schematic picture of the in situ observation of SA 
by YM. White and black circles denote gold and copper 
atoms, respectively, (a) Before the onset of SA, individ- 
ual copper atoms are deposited on a gold cluster, (b) 
After SA is completed, copper atoms dissolve into a gold 
cluster to form a homogeneously mixed alloy cluster. 

They concluded that the remarkable features of SA 
phenomenon can be summarized as follows: 
(l)The diffusion rate of copper atoms in clusters is 
approximately 9 orders of magnitude faster than that 



The experimental results mentioned above suggest the 
presence of an unexpectedly fast diffusion process, which 
is controlled by negative heat of solution and tempera- 
ture, in a small sized cluster. 

It is worth while to state what is unusual in SA. It 
is quite natural for Au and Cu atoms to mix and and 
change into an alloyed state, which is entropically prefer- 
able from the viewpoint of equilibrium statistical me- 
chanics. In addition, negative heat of solution ensures 
that the alloyed state is enthalpicaly preferable. In fact 
near icosahedral-like ground states, mixed structures are 
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preferred over the segregated ones for bimetallic clusters 
with 13—14 atoms |1. Thus it seems to be no sur- 
prise to find out a spontaneously alloyed state in a binary 
cluster. The really unusual point of SA is that the al- 
loying completes at least within the time scale of second 
even at room temperature. How is this diffusion real- 
ized? What we would like to elucidate is the atomistic 
mechanism of such a fast diffusion. Before pursuing this 
question, it is helpful to confirm how fast the diffusion 
is. Supposing that a diffusion process of impurity atoms 
obeys Arrhenius law, YM's experimental result implies 
that the activation energy is effectively lowered, at least, 
by 40—50% in Au clusters compared with bulk crystalline 
Au. ft is important to note that the diffusion coefficient 
of atoms on a clean surface is many orders of magnitude 
faster than that in bulk. In fact the activation energy of 
the diffusion constant on a clean surface is about 35% of 
that in the bulk ||. That is, the diffusion coefficient of 
Cu atoms in a nano-sized cluster is almost comparable to 
that on a surface. Atoms are easy to move on a cluster 
surface, because the surface is populated by point defects. 
Considering that the fraction of surface atoms in a cluster 
becomes significant as the cluster size decreases, it is clear 
that the rapid surface diffusion is relevant for the rapid 
SA phenomenon. On the other hand, it is important to 
note that SA is a manifestation of atomic movement in 
the radial direction of the cluster. How is the fast diffu- 
sion parallel to the surface of cluster related to the rapid 
diffusion in the radial direction? A remarkable difference 
between a bulk surface and a cluster surface is that the 
latter has a non-vanishing curvature, which makes the 
surface more deformable. An easily deformable surface 
may influence the manner of atomic diffusion along the 
surface and may drastically accelerate the radial diffusion 
rate ||lo| . How is such a rapid diffusion enhanced or sup- 
pressed by the key factors such as magnitude of negative 
heat of solution, cluster size and temperature rise? 

In our previous work, we reported some preliminary 
results which included the numerical simulation of SA 
by using a 2-dimensional(2D) binary Morse model [fl3|| . 
We pointed out that some features of the experimental 
results summarized by (l)-(5) may be reproduced. How- 
ever, the results presented there were not complete. The 
time scale of simulation was not long enough to repro- 
duce the whole SA process in larger-sized clusters and 
the tempareture of simulation was not lower enough than 
the melting temperature. In particular, the dependence 
of the alloying process on the cluster size, which is the 
heart of YM's experiments, was not made clear. The aim 
of the present paper is to show systematic results of a fur- 
ther extensive and comprehensive numerical simulation 
and demonstrate how our simple-minded Morse model of 
clusters reproduces the essential feature of YM's exper- 
iment. A particular emphasis is put on the fluctuation 
and rearrangement of cluster atoms, which contain use- 
ful information for elucidating the atomistic mechanism 



of the SA process. 

III. A NUMERICAL APPROACH: HOW TO 
PREPARE A MODEL WITH APPROPRIATE 
INITIAL CONDITIONS 

A. A model cluster 

There arc several proposed empirical potentials which 
successfully mimic the equilibrium properties of bulk 
metals such as a lattice constant, bulk modulus, elastic 
constants and sublimation energy. Among them the so- 
called Embedded Atom Method(EAM) is a well-known 
model for alloy systems However, it has many pa- 
rameters which should be well-adjusted to yield plausible 
values for equilibrium properties of bulk metals. More- 
over, transferability or applicability of these potentials to 
a cluster system is still unknown [H . Because a simpler 
model is better to get physical insight into the mechanism 
of SA, we employ a Morse model which has can explain 
qualitative aspects of the experimental results jl3];[l4} . 

More specifically, we choose the Morse model for the 
two reasons: 

(1) Unlike EAM, the heat of solution, which is the key 
parameter of our simulation, can easily be controlled by 
a single parameter. EAM has many parameters which is 
influential in changing heat of solution. 

(2) By using pairwise potential we can considerably re- 
duce the simulation time. 

Furthermore, in the present work we use a 2D Morse clus- 
ter rather than the realistic 3D cluster. The reasons why 
we examine the 2D model are twofold. First, the compu- 
tation time for the 2D model is much shorter than that 
for the 3D model. In realistic 3D simulations correspond- 
ing to the experimental condition, time evolution of more 
than 1000 atoms should be traced for longer than a few 
microsecond. This is still beyond recent computational 
capability. Secondly, visualization of individual atomic 
motion in the SA process can more easily be done with 
2D model than that with 3D model. 

We take 2D Morse potential Hamiltonian 

N 

i— 1 i<j 

and 

V k i(r) = e kl {e- 2Mr - r ^ - 2e~^ r - r ^}, (2) 

where k and I specify the two species of atoms, say host 
and guest. Host and guest atoms are denoted by A and 
B, respectively. A cluster is formed by Na host atoms 
and Nb guest(or solute) atoms, where total number of 
atoms is TV = Na + Nb- For simplicity we choose Paa = 
Pbb = Pab = 1.3588IA- 1 ], £ = £.4.4 = e BB = 0.3429[eV] 
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and r AA — r c BB — r c AF , — 2.866L4]. Those values are 
suitable for copper Jig ]. The only free parameter is 
a = cab/caa- Because the heat of solution AH is given 
by AH = z(l—a)e where z is a coordination number, our 
choice for a (a = 1.1), provides a negative heat of solution 
for the binary system Jl6| . For a realistic binary system 
the relations, eaa = (bb = £ab and r c A 



AA 



BB 



AB 



do not hold. Although we may oversimplify the model 
system, we believe our model is suitable to investigate 
how the effect of the heat of solution controls SA. 

In our model we neglect the presence of supporting film 
which may be important as a heat reservoir. However, 
the released binding energy is transfered to a substrate 
very slowly, since the coupling between a substrate and 
a cluster is considerably weak. Weak coupling between 
a substrate and a cluster is due to the frequency mis- 
match between atoms in a substrate and a cluster(See 
appendix) . 



B. Plausible initial conditions 

It is very hard to simulate the realistic condition of 
YM's original experiment. The solute atoms which are 
successively deposited onto the host cluster form new 
bondings and release the bonding energy as excess ki- 
netic energy. This process makes the cluster so hot that 
some cluster atoms suddenly evaporate, because the to- 
tal energy is conserved in our simulation. We, therefore, 
employ an initial condition in which the solute atoms 
are bounded stably with the host cluster at its surface 
as shown in Fig. 2 in order to remove other initial con- 
ditions in which surface atoms are suddenly evaporated. 
In addition, as will be discussed later in detail, it is im- 
portant to note that the temperature, which is identified 
with the average kinetic energy of the cluster, is one of 
the key parameters which controls SA. 
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Fig. 2 A typical initial atomic configuration for isoen- 
ergetic MD simulation. White and black circles denote 
host(A) and guest(B) atoms, respectively. A cluster con- 
sists of 47 host atoms and 20 guest atoms, (Na, Nb) = 
(47,20). (The radius of each circle is given by an arbi- 
trary unit.) 



Thus it is desirable to keep the initial temperature of 
the cluster as stationary as possible. Arbitrarily chosen 
initial configurations are unstable in the sense that the 
surface atoms do not form saturated bondings with the 
inner atoms and the configuration is not located in a suf- 
ficiently deep local minimum of the potential energy. As 
a result, in time evolution of MD run, formation of new 
bonds between atoms leads to a rapid and considerable 
temperature rise. In such cases the initial temperature is 
not well controlled. Indeed, we have observed, in prelim- 
inary simulations, that atoms attached to the surface of 
a cluster are in general promptly absorbed into the first 
layer of a host cluster and release the kinetic energy in 
an uncontrollable manner. These are the practical rea- 
sons why we choose a stable configuration as is shown 
in Fig. 2 with which the kinetic energy keeps a station- 
ary value over a sufficiently long time scale in which the 
initial temperature is well defined. 

Selection of a stable initial configuration is also jus- 
tified from the experimental viewpoint. Using nm-sized 
amorphous Sb(a-Sb) clusters, each of which was attached 
to a Au cluster, YM also confirmed the presence of SA 
|l8| . In their experiments these binary clusters were 
gradually heated from 96K to 290A". Then, dissolution 
of Au into a-Sb clusters set in around 200K and Sb-Au 
alloyed clusters were produced in the time scale much less 
than lsec. Unlike the original experiment of (Au — Cu) 
cluster by YM, the two sorts of clusters bonded firmly 
with each other at their interface and could be identified 
with our stable initial configuration. 

Stable initial configuration for MD is prepared by the 
simulated annealing method: prepare a homogeneous 
cluster composed of the A-atoms alone and start its time 
evolution, then the system wanders over various local 
minima of the potential energy. At the step when the 
configuration falls into a sufficiently deep local minima, 
we choose it as the initial condition and start the sim- 
ulation of SA after assigning an appropriate number of 
atoms at the surface of the cluster to the solute atoms 
i.e., B-atoms. After allocating initial random velocity to 
each atom, the trajectories were computed by the veloc- 
ity form of the Verlet algorithm, where time step used 



was 5 x 10 °sec [Jit 
run in every case. 



Time evolution of 10 sec was 



IV. TOOLS FOR QUANTIFICATION: HOW TO 
OBSERVE ATOMIC FLUCTUATION AND 
REARRANGEMENT DURING SA 



In this section we give some examples of numerical re- 
sults obtained from microcanonical MD simulation. Our 
main purpose here is to introduce some useful quantities 
and to outline how SA proceeds by them. 
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A. Useful tools 

Prior to discuss numerical results of size, temperature 
and a dependency, we introduce some quantities char- 
acterizing the atomic motion in a cluster. At first, we 
define a quantity measuring to what extent two different 
types of atoms are mixed in a cluster. To do this, we 
introduce the number of neighboring host atoms per a 
solute atoms, say nsit). is defined as, 



1 



N B 



i—1 



(3) 



where N A m is the number of A-atoms which occupy 
the nearest neighbor sites around the i-th B-atom at 
time t. In case that two types of atoms are homoge- 
neously mixed, a simple mean field consideration yields 
that n,B{t) should be n 1 ^ — z(l — r), where z and r de- 
note coordination number and fraction of solute atoms 
respect to total number of atoms in a cluster. 

During the SA process each atom in a cluster vibrates 
near a stable position and sometimes jumps to another 
neighboring position. The accumulation of the latter pro- 
cess results in the mixing of the solute atoms into the host 
atoms. Accordingly the dynamics of alloying process has 
at least three different time scales: the shortest one is due 
to the rapid fluctuations around the site, which is com- 
parable to the inverse of Debye frequency ( ~ 10 _1 [ps]). 
The second one, which is much longer than the first one, 
is the time scale of rearrangement of their neighboring 
atoms (~ 10 — 100 [ps]). The longest one characterizes 
the time scale of the alloying process which is the relax- 
ation time from the initial nonequilibrium configuration 
to an equilibrium one, which is shorter than 10msec in 
YM's experiments. Since the third time scale can be ob- 
served in terms of the variation of ns(t), we introduce 
alternative quantities characterizing atomic motion ob- 
served during the shorter two time scales. To quantify 
fluctuating and rearranging properties of atoms, we in- 
troduce the atomic Lindemann index 5(i) and the fre- 
quency of recombination of the neighboring atoms. The 
rapid vibration of atoms around each site is manifested 
in the so-called Lindemann index, which is expressed by 
the root mean square(rms) deviation of atomic separation 
between neighboring atoms ||. We define the nearest 
neighbor Lindemann index(NNL) for individual atoms, 
5i(t), as follows; 



Si(t) 



1 



E 

jGn.n. 
of i — th atom 



(Ri. 



(R 



ij/t 



(4) 



where Rfj denotes the distance between i — th and j — th 
atoms, (N^')t is time-averaged number of the nearest 
neighbor atoms of i — th atom. Note that (T) is time 
average of the arbitrary quantity T, given by, 



to 



t+t a 



T{t)cLt. 



(5) 



The averaging time t av is fixed to be 2 ns. On the other 
hand, the frequency of recombination of the neighboring 
atoms is estimated by the distance index Q. Distance 
index is derived from a adjacency matrix, say M, which 



is N x N symmetric matrix whose elements 



= 1 for 



I r ij |< r c and zero otherwise, where r c = l.?>Ar c AA . Dis 
tance index di(t) of the i-th atom is, then defined as 



di(t) = 



\ 



N 



I My(t + At)-My(t) I 1 



(6) 



Supposing that atoms move from one site to another site 
frequently in a cluster, then the occurence of atomic rear- 
rangements should be manifested by the variation of the 
distance index di(t). Time interval, At, must be short 
enough to resolve the single event of atomic rearrange- 
ment. In our numerical analysis we set At = lOps. 



B. Some examples 



Before illustrating systematic results of our simula- 
tion, we show a typical example of time evolution of 
the alloying process observed in our numerical simu- 
lation for the A47B20 cluster (denoting the cluster of 
(N A ,N B ) = (47,20)) and a = 1.1. In Fig.3 we show 
the time evolution of 71b(£), i.e., the number of neighbor- 
ing atoms of different species per a solute atom, nsit) 
increases monotonically from the initial value tib(0) ~ 2, 
which means that the solute B-atom initially on the sur- 
face of the cluster forms two bonds with the host A-atoms 
in the inner shells. Finally ns(t) reaches to the level 
of ng corresponding to the homogeneously mixed state, 
which is indicated by the dotted line. 




500 

Time [ns] 
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Fig. 3 A typical time evolution of n B (t) for isoenergetic 
MD sumulation. The final configuration of atoms in 
the cluster is also inserted. Initial temperature is about 
510[K].(See Fig.4.) 

It is obvious that a homogeneously mixed configura- 
tion is almost achieved until 800 [ns]. The present ex- 
ample corresponds to the case where the initial temper- 
ature is the lowest among our data exhibiting SA within 
1000 [ns]. The spontaneous increase of the Us(t) implies 
that the system evolves so as to decrease the potential en- 
ergy of the system. Indeed, the increment ,n B (t) — n B (0), 
means that the variation of potential energy per one 
B-atom is (n B (t) — n B (0))z(l — a). The time aver- 
aged variation of the total potential energy, then, is 
(AU(t)) t = N B (n B (t) - n B (0))(l - a), if we take into 
account the contribution from the nearest neighboring 
atoms. The decrease in the potential energy is converted 
into the increase in the kinetic energy. In our simula- 
tion we define the kinetic temperature by T = k B (2N~z) > 
where Ek is total kinetic energy of the system and k B is 
Boltzmann constant. Note that we eliminate the contri- 
butions from translational and angular degrees of free- 
dom, because we select initial conditions with vanishing 
translational and angular momentum. Variation of the 
kinetic temperature of the cluster is shown in Fig.4. 




500 



Time [ns] 



Fig.4 A typical time evolution of kinetic temperature 
for isoenergetic MD simulation of A47B20 binary Morse 
cluster taken from the same data as Fig. 3. Initial tem- 
perature is estimated as 510 [A"]. Averaging of kinetic 
temperature was done in every 2[ns]. Notice that re- 
sulting temperature does not exceed estimated melting 
temperature(~ 670 [AT]) (See also Table. I). 

As is expected, it is clear that the variation of n B (t) 
is strongly correlated to the temperature rise. The fi- 
nal increase in the temperature, say AT, is less than 



AT max = - kB{2 2 N -3) ( A U{t = oo)) tav ~ 200[K], and 
the half of AT max contributes to the actual tempera- 
ture rise according to the virial theorem. Indeed the 
virial theorem predicts AT ~ 100A, which is consistent 
with Fig.4. Although the temperature is increased up 
to about 610 [K] when the alloying process is completed, 
the temperature at t = 800 [ns] is still sufficiently below 
the melting temperature of the cluster. The typical melt- 
ing temperatures measured for various sizes of clusters is 
listed in Table I. 

Although one cannot observe a sharp solid-liquid tran- 
sition in a small system in a strict sense, we can prac- 
tically locate the melting point by an abrupt jump in 
caloric curve and Lindemann index |^,^3|. As is well- 
known experimentally and numerically, the melting point 
is reduced as size of cluster decreases Q . It is evident 
that Table I also indicates the same trend. 

Table. I The relation between cluster size and the melt- 
ing temperature for our 2D Morse model. 



Number of atoms 


Melting tempetature [K] 


32 


580 


67 


670 


80 


710 


117 


750 



Since the melting temperature of a cluster A47B20 is 
about 670 [K], dynamics of SA process in Fig.4 provides 
an evidence manifesting that the alloying process com- 
pletes in the solid phase without the melting of the whole 
cluster. 

In Fig. 5 we show the snapshots of the atomic config- 
urations corresponding to SA process depicted in Fig. 3 
and 4. 

•mo*o m j-^to*.-, "•••£o. *°o»°o. 
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580ns 610ns 755ns 880ns 

Fig. 5 Snapshots of atomic configuration of A47B20 in 
time evolution. White and black circles denote host and 
guest atoms, respectively. Shaded circles are host atoms 
which initially forms core of a cluster. 
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The B-atoms initially deposited on the surface wan- 
ders actively along the surface of the cluster. However, as 
clearly seen in Fig. 5, they all stay on the surface and the 
movement in the radial direction of the cluster is almost 
quiescent over a significantly long time (t < 300[nsec]). 
The length of the quiescent period depends very sensi- 
tively on the initial configuration. After that the radial 
movement begins to be activated. Some atoms suddenly 
enter into the inner shell, which is reflected in an abrupt 
rise of ns{t) in Fig. 3. Then the B-atoms enter further 
into inner layers and ns{t) increases stationarily. During 
the stationary stage, the B-atoms aggregate in one side 
of the cluster as is typically seen in the atomic config- 
urations at t = 520 — 580 [nsec]. The atoms initially 
forming the core of the cluster, which are indicated by 
shaded atoms in snapshots, arc pushed out in a group 
and breaks up when they reach the surface. In this way 
the outer B-atoms exchange their positions with the in- 
ner A-atoms. All the alloying processes observed in our 
microcanonical simulation completes according to the 
similar scenario. 



V. DEPENDENCE OF THE ALLOYING 
PROCESS ON KEY PARAMETERS 

According to the conclusions by YM, SA phenomena 
are dominated by three key parameters; magnitude of 
negative heat of solution, temperature of the system and 
the size of a cluster. In this section we show system- 
atic numerical results exhibiting that these parameters 
are similarly important for the onset of SA reproduced 
in our microcanonical simulation. 



A. A role of heat of solution as a driving force of SA 



mixed to some extent. However, it does never reach the 
value, rig, within the simulation time(800[ns]). The dif- 
ference in the variation of ns{t) is a direct evidence in- 
dicating that the SA process is dominated by heat of 
solution. In other words, time to complete SA becomes 
shorter as initial temperature is getting higher. A sys- 
tematic analysis to clarify the relationship between the 
alloying speed and initial temperature is pursued in Sec 
V.B. 
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In Sec. IV we illustrated a typical dynamical behav- 
ior of SA which is driven by negative heat of solution. 
In order to probe the dependence of alloying process on 
heat of solution we compare three different cases, i.e., 
a = 0.9(positive heat of solution), a = 1.0(vanishing 
heat of solution) and a = 1.1 (negative heat of solution). 

In Fig. 6 (a) (b), we show the time evolution of ns(i) 
for two values of initial temperature, 620[if] and 510[K], 
respectively. In both cases we also display the time evo- 
lution of ns(t) for a = 0.9, a = 1.0, and a = 1.1. 

Some significant discrepancies in the variation of ns{t) 
among these three cases become obvious in Fig. 6 (a)and 
(b). It is clear that the value of ns{t) for a = 1.1 shows 
a rapid increasing trend which is a signature of a faster 
alloying process. Conversely, an absence of a mixing be- 
tween guest and host atoms is manifested in a slower de- 
crease and saturation of ns(t) for a = 0.9. For a = 1.0, 
rig (t) increases very slowly, because A- and B-atoms are 



Fig. 6 A typical example of time evolution of ns{t), 
which is defined for host atoms (A-atoms) and guest 
atoms (B-atoms), obtained by isoenergetic MD runs with 
high and low initial temperature. ((a)610[K](high) and 
(b)510[K](low)). (c)The radial diffusion of atoms in a 
cluster is evaluated in terms of the time evolution of the 
distance between B-atoms and the center of mass of the 
cluster for the case(a). The final configurations of a clus- 
ter for a = 1.1 and a — 1.0 are also inserted. An ag- 
gregation of B-atoms in a cluster is also indicated by the 
circle of the inserted figure for a = 1.0. 

In addition to the main features mentioned above, it 
is important to note the following numerical results. 
(l)For the relatively high temperature, the mixing be- 
havior takes place even for the atomic species with null 
heat of solution (a = 1.0). As shown in Fig. 6(c), most 
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of the B-atoms diffuse into the cluster even for the case 
of a = 1.0. That is, the average distance between the B- 
atoms and the center of mass reaches approximately 8.8, 
which is very close to the value of the final configuration 
of a — 1.1, where a complete mixing is attained. 
(2) On the contrary, ns(t) for a — 1.0 is always smaller 
than that for a = 1.1. 

These seemingly contradictory facts are understand- 
able by comparing the final atomic configurations for 
a = 1.0 with that for a = 1.1. As shown in Fig.6.(c), the 
A-atoms and the B-atoms are well separated for a = 1.1, 
while the B-atoms tend to aggregate for a = 1.0. Aggre- 
gated configurations of the A- and the B-atoms are en- 
ergetically neither favorable nor unfavorable for a = 1.0, 
although such an aggregation of the B-atoms are ener- 
getically unfavorable for a > 1.0. Indeed, as indicated 
by the circle in the inserted figure of Fig. 6(c), an aggre- 
gation of the B-atoms is easily verified. In short, in case 
of a = 1.1, the A- and the B-atoms mix so as to decrease 
the potential energy and it works as the driving force 
of SA. Consequently, the role of negative heat of solu- 
tion is twofold, as far as our microcanonical simulation 
is concerned. First, it results in a driving force which 
promotes the alloying of two species of atoms. Secondly, 
the decrement of potential energy in the alloying process 
is converted into the kinetic energy, which heats up the 
cluster and accelerates SA. The latter effect is persued 
again in the next section. However, these two facts in- 
dicate that the diffusion of the B-atoms into the cluster 
from the surface is by no means prohibited even in case 
of a — 1.0, where the potential energy gain due to the 
mixing is zero. It should be noted that such a diffusion 
process with is still much quicker than the diffusion into 
the bulk media. We numerically confirmed that no sig- 
nificant diffusion occurs in the bulk 2D medium within 
the time scale of Fig. 6 at the same temperature. 

In summary, heat of solution is a key parameter of SA 
in the sense that positive heat of solution prohibits SA, 
while negative heat of solution remarkably accelerates 
SA. However, a rapid mixing of the two species of atoms 
occurs even in case of null heat of solution. This fact 
demonstrates that the rapid diffusion process is a generic 
feature of microclusters. The present result is consistent 
with YM's experiments. Indeed, YM reported that SA 
occurs even in the combinations of atomic species with 
very small magnitude(almost null) of positive heat of so- 
lution, when the size of cluster is sufficiently small pjj , 

B. Temperature dependence and size dependence 



canonical dynamics, we regard the initial kinetic temper- 
ature as the parameter characterizing the temperature 
of the system. The temperature dependence of the alloy- 
ing time is examined for the four sizes of clusters A23B9, 
A4JB20, ^56^24, and Ai4oBqq. For each of them, we 
prepare initial conditions corresponding to various val- 
ues of initial kinetic temperature. The values of initial 
kinetic temperature cover a wide range from well below 
to just below the melting temperature. The value of 
the paramter a is fixed (a = 1.1) for all sizes of clus- 
ters. Note that the ratio r = N . B , r is chosen to be the 
same value, r = 0.3, for the four sizes of clusters. This 
is because the bonding number of homogeneous mixing 
rig should be fixed to be a common value irrespective of 
the cluster size. Moreover, the possible temperature rise 
AT = 2^[ n B ~ riB(0)]r(a — 1) in the SA process are 
controlled to be common in the four sizes. 

The alloying time should be defined as the time re- 
quired to attain the homogeneous mixing, which cor- 
responds to the time when nsit) reaches the value for 
homogeneous mixing i.e.,n|| = 4.2. However, it takes 
extremely long time to complete SA especially for larger 
clusters, and it is difficult to get reliable data. We, there- 
fore, define the alloying time T a u oy as the time when 
tib (i) reaches 3.0, which is almost the average of the ini- 
tial value of tib and rig. The semi- log plot of the inverse 
initial temperature versus the alloying time is depicted 
in Fig. 7 for the four clusters. 




6.0015 ' 0.002 6.0025 
1/T [1/K] 



In YM's experiments, the time needed for the SA pro- 
cess to complete depends sensitively on the temperature. 
Some systematic results for the initial temperature de- 
pendency of SA is presented here. Because kinetic tem- 
perature is not constant during time evolution in a micro- 



8 



Fig. 7 The dependence of the alloying time T a u oy 
on the initial temperature To for 2D Morse binary 
cluster A23-B9 (black square), A47B20 (white square), 
^56^24 (white circle), and Ai4 B 60 (black circle). 

Because SA is an outcome of the diffusion of solute 
atoms into the cluster, it is reasonable to expect that the 
alloying rate also obeys the Arrhenius-like law of the dif- 
fusion coefficient. Since the time evolution pattern ns(t) 
largely fluctuates from sample to sample, the alloying 
time defined above is also accompanied with consider- 
able sample-dependent fluctuations. However, as shown 
in Fig. 7 the semi- log plots for these samples are almost 
on lines. Aftrer all, the dependence of the alloying time 
on the initial temperature obeys an Arrhenius-like law: 

Talloy = To exp(— 7— ), (7) 

Jo 

where T a u oy can be interpreted as the 'activation energy' 
divided by Boltzmann constant and tq is the inverse of 
the frequency factor. It is evident that the activation 
energy, which is the slope of each line, increases with 
the size of the cluster. The present result is the mani- 
festation of the size effect which was observed in YM's 
experiments and is regarded as the characteristic effect 
peculiar to clusters. The activation energy T a ii oy and the 
inverse of the frequency factor To evaluated for the four 
sizes of clusters are summarized in Tabic II. 

Table II Numerically estimated To and T a iio y for various 
sizes of binary clusters, i.e. Ai4 B &0 , A 5 qB 2 4 1 ^47^20 and 
A23B9. 



Cluster Size 


to [psec] 


Talloy [K] 


A\4oBqo 


0.02 


8600 


^56^24 


2.7 


5400 


A47B20 


0.91 


6000 


A23BQ 


2.7 


3000 



By an extrapolation with these values we roughly es- 
timate the alloying time at room temperature. The esti- 
mated alloying times are 5.4 x 10~ 10 [sec], 7.3 x 1.0 _4 [sec], 
4.4 x 1.0~ 4 [sec] and 1.3 [sec] for A 23 B 9 , A 4r B 2 o, A 56 B 2 4 
and Ai4 B 60 at 300K, respectively. The resulting values 
are short enough to be consistent with YM's experimen- 
tal observation. 

In addition, it is interesting to note the fact that ac- 
tivation energy T a iio y has similar values for A 56 B 2 4 and 
A47B20, although the size of the former is larger than 
that of the latter. The apparent contradiction is under- 
standable, if we assume that these clusters consist of the 
same number of shells. (The precise definition of the shell 
is given in Sec. VI. A.) One can easily verify that the clus- 
ters, A23-B9, A47B20, A 5( jB 2 4 and A^Bqo are composed 
of 2,4,4, and 7 shells, respectively. Thus, the cluster size 
represented by the number of shells is a relevant quantity 
to determine the activation energy given in Tabic. II. 



VI. AN ACTIVITY OF A CLUSTER SURFACE 
AND THE MECHANISM OF SA 



A. A shell dependent activity of atoms 



In this section we probe how actively individual atoms 
composing the cluste move during SA. To evaluate ac- 
tivity of atomic motion in a cluster provides important 
clues to elucidate the mechanism of SA process. In par- 
ticular, we pay our attention to how dynamical activities 
represented by fluctuation and rearrangement of atoms 
depends on the distance from the center of cluster. To 
do this, it is convenient to divide the cluster into shells. 
According to the distance of a target atom from the cen- 
ter atom, which is defined as the atom closest to the 
center of mass of the cluster, we allocate the shell in- 
dex to each atom. Since a single cluster has a hexagonal 
structure, it is possible to introduce magic number where 
a cluster forms a geometrically packed configuration. For 
instance, a cluster which contains 7, 19, and 37 atoms are 
magic number clusters which consist of 1, 2, and 3 closed 
shells, respectively. A cluster containing 67 atoms, which 
is shown in Fig. 3, is divided into 4 shells, say m = 4. The 
shell index number m is assigned in order of the distance 
from the center of mass. The center of mass atom is 
allocated to the zero-th shell. In Fig. 8 the frequency dis- 
tribution of the distance of atoms from the center of mass 
is depicted. This is a typical example obtained from a 
single isoenergetic MD run of a A47B20 cluster. 




Distance from cm [A] 



Fig. 8 A frequency distribution of a distance between the 
center atom and other atoms in a 2D Morse binary clus- 
ter A47B20 ■ The dotted line denotes critical distances 
used to divide a cluster into shells. 

Since the distribution has distinct minima in Fig. 8, it is 
possible to divide a cluster into shells without any redun- 
dancies. In the followings we argue about the fluctuation 
and the rearrangement of atoms in an individual shell. 
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A fluctuation property of the separation between 
nearby two atoms has been used as an indicator to lo- 
cate the melting point ||. In fact, the location of the 
melting point determined by rms bond length fluctu- 
ation(Lindemann index) almost coincides with the one 
given by the caloric curve |^|,^2)- In order to get fur- 
ther detailed information on atomic fluctuation we use 
the nearest neighbor Lindemann index(NNL) defined for 
each shell. The NNL index of the k — th shell, Afc(t), is 
defined by averaging the NNL index for individual atoms, 
say 5i(t), belonging to the same shell; 



A fc (t) 



where {N<p ) is time average of the total number of atoms 
contained in the k — th shell over t av . 

One can expect that some shell-dependent dynamical 
activities are captured by A/-(i). In Fig. 9 we depict the 
time evolution of Afc(t) for a single run of SA process. 
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Fig. 9 (a)Time evolution of Lindemann index for the 1st, 
2nd, 3rd and 4th shell obtained from data of 2D Morse bi- 
nary cluster A^jB2q- (b)Time evolution of the cumulated 
distance index for (k = 1 — 4). Initial temperature 
is about 510[K].(See Fig.4.) 

The bond fluctuation is much enhanced in the shells 
near the surface of a cluster even if its temperature is 
substantially below the melting point, which is about 
680 [K]. It exceeds the Lindemann criterion for melting 
i.e., Afe ~ 0.1, denoted by the broken line in Fig. 9. The 
time evolution of A4 in Fig. 9 (a) strongly suggests that 
the cluster surface is in a melting state. Indeed the sur- 
face melting is observed in Pb cluster below melting tem- 
perature [ fL0| . Judging from the fact that the Lindemann 
index of the inner shells are less than 0.1, the inner shells 
of the cluster is in a solid phase in the sense of Linde- 
mann's criterion. 

The dependence of the Lindemann index for the 3rd 
and the 4th shell upon the initial temperature is shown 
in Fig.10. 
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Fig.10 The relationship between the initial tempera- 
ture and the averaged atomic fluctuation for the 3rd 
and the 4th shell in a 2D Morse binary cluster A^jB2q- 
The averaged atomic fluctuation is given by A& = 
TTi — fo' 1 "° y dtA k (t). The black and white circles denote 
the Lindemann index for the 3rd and the 4th shell, re- 
spectively. The surface atoms on the 4th shell begin to 
melt near the region denoted by the arrow. 

While the outer most shell begins to melt about 510[K], 
SA completes within 800 [ns] at the same temperature as 
shown in Fig. 3 and 4. Consequently, it is plausible to say 
that the presence of the active surface atoms is necessary 
condition to attain SA within Ifisec. 

The Lindemann index measures the fluctuation occur- 
ring on the time scale of atomic vibration, which is the 
shortest time scale of the system. An enhancement of 
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fluctuation on the surface also increases the chance for 
the surface atoms to jump to another vacant sites on the 
surface. Such a rearrangement process is a rare event 
which occurs in time scale much longer than the vibra- 
tion. On the surface, however, the atomic rearrangement 
is also enhanced, because the surface is filled with de- 
fects or vacancies to which surface atoms are allowed to 
jump. Ease of jump on the surface enhances the diffu- 
sion constant along the surface. For a demonstration of a 
pecular motion of surface atoms, it is interesting to show 
how the activity of atomic rearranging process depends 
on the shells. In order to quantify it in each shell we 
define the cumulated distance index for the fc— th shell, 



C%\t)= ( drD k ( T ), 
Jo 



(9) 



where 



W) = ^y £ *(*)■ (10) 



£k — th shell 

where di(t) is the distance index introduced by eq.(6), 
and Cp\t) is the accumulated number of rearranging 
events which occur near the atoms belonging to fc— th 
shell. As easily verified in Fig. 9(b), atomic rearrange- 
ment is much more frequent in the surface shell, and al- 
most all the rearranging events occur in the shells on or 
close to the surface. 



^4-56^24, we obtain Tr ~ 3700 [K] for the surface shell. 
On the other hand, as shown in Table. II, the activation 
energy of alloying, T a u oy , is about 5400 [K]. These two 
activation energies are significantly different, 



' alloy 

~Tr~ 



1.5. 



(13) 



It is possible to interpret Tr as the activation energy of 
rearranging motion parallel to the surface, where T a u oy 
measures the activation energy for the atoms to diffuse 
along the direction perpendicular to the surface. Let Pu 
and Pj_ be the probabilities for surface atoms to jump 
to another site parallel to and perpendicular to the sur- 
face, respectively. Supposing that the radial hopping is 
induced by the parallel rearrangement, the relation be- 
tween these two probabilities is expressed as, 



_£ = e -T R p/T 



(14) 



where Trp — {T a iioy — Tr) > 0. A plausible interpre- 
tation of Eqs.(13) and (14) is that the diffusion in the 
radial direction of cluster is the outcome of a surface 
rearrangement followed by a certain activation process 
characterized by the barrier height T RP , which is roughly 
evaluated as \T R . The present interpretation is also sup- 
ported by a direct observation of trajectories of atoms 
during the SA process. In Fig. 11 the trajectories of all 
solute atoms (B-atoms) are shown for every 250[ns]. 



B. Surface melting and the rapid radial diffusion 



It is possible to demostrate that the surface activity 
is responsible for the SA process in an alternative way. 
Suppose that we could suppress the activity of surface 
atoms, for example, by embedding the alloying cluster in 
a bulk medium, then the rapid alloying process would be 
inhibited because the diffusion in the radial direction of 
the cluster is nothing more than the diffusion in a bulk 
medium. The active rearrangement of the surface atoms 
is a necessary condition for the rapid alloying to be re- 
alized. The frequency of the rearrangements per unit 
time is represented by Dk(t), which increases as the time 
elapses due to the temperature rise. Its average over a 
single alloying process, say D k , obeys Arrhenius-type law 
respect to initial temperature Tq: 

T 

D k oc exp[— (11) 

where 

D k = / dt'D k (t'). (12) 

Talloy JO 

This nice property enables us to introduce the activation 
energy Tr of the atomic rearrangement. For the cluster 




Fig. 11 Trajectories of solute atoms (B-atoms) during SA 
in a cluster A/^iB2o{ot = 1.1). Initial configuration and 
initial temperature are the same as those in Fig. 5. and 
Fig. 4, respectively. 

During the initial stage of the process (<~ 260 [ns]) the 
B-atoms glide over the surface of a cluster. In the next 
stage the atomic motion begins to contain the compo- 
nent perpendicular to the surface and a diffusive motion 
into the cluster is activated^ 510[ns]). The trajecto- 
ries of B-atoms spread over the whole cluster in the final 
stage. In short, atomic gliding motion along the surface 
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plays a role of a trigger to push solute atoms, which was 
initially located on the surface, into the inside of a clus- 
ter. Successive snapshots in Fig. 11 demonstrate that the 
diffusion along the surface induces the diffusion to the 
inside of cluster. 

In addition, the diffusion into the solid cluster, which 
is responsible for the SA process, is quite different from 
the diffusion into the bulk solid in the following aspect. 
A usual diffusion mechanism of atoms into the bulk is 
attributed to a hopping process which is mediated by 



defects or vacancies in the solid 23 1. In contrast, a for- 
mation of defects in clusters is an extremely rare and its 
lifetime is very short, because defects are immediately 
pushed out to the surface. Then the resulting hopping 
probability in the radial direction mediated by the de- 
fects is negligibly small. Moreover, in our simulations, 
we could find no evidence suggesting that the vacancies 
inside the cluster play any role in the diffusion process 
in the radial direction of the cluster. Thus we conclude 
that the diffusion into the solid cluster, which is relevant 
for the SA process, is quite different from the diffusion 
process in the bulk solid. To elucidate the essential mech- 
anism of the rapid diffusion into the solid cluster will be 
discussed elsewhere in detail ^4j. In the present work, 
we only claim that, as discussed above, the frequent rear- 
rangement of atoms on the surface of the cluster initiates 
the rapid diffusion in the direction perpendicular to the 
surface. 



VII. CONCLUSION 

In the present paper we examined the presence and the 
absence of SA in a 2D binary microcluster in terms of an 
isoenergetic MD of a simple Morse model under different 
conditions of initial temperatures, sizes, and heat of solu- 
tions. One important advantage of our model is that the 
maginitude of negative heat of solution, which has been 
supposed to be he key parameter of SA in YM's experi- 
ment, can be well controlled by a single parameter. We 
confirmed that the heat of solution is similarly a driving 
force to form homogeneously mixed binary cluster in a 
very short time scale less than l[/is] in our model. Our 
main conclusions consist of the following three results. 

(1) By changing the initial temperature of the system, 
we found that SA occurs sufficiently below the melt- 
ing temperature. The time required to complete SA 
becomes longer exponentially as the initial temperature 
To decreases. That is, it obeys an Arrhenius-like law, 

Talloy OC exp[-^yr^]. 

(2) We investigated the cluster size dependence of the 
alloying time. As a result, we found observed that larger 
size cluster spend longer time to achieve SA. More pre- 
cisely, the activation energy T a u oy becomes larger with 
increase in the cluster size. This result makes it clear 
that the quicker alloying surely occurs in smaller sized 



clusters. By extrapolating the Arrhenius plot, we veri- 
fied that the alloying time is much less than sec order at 
room temperature for a sufficiently small cluster. These 
numerical results qualitatively coincide with the experi- 
mental observation by YM. 

(3) By introducing quantites to probe fluctuating and 
rearranging properties of atoms composing a cluster, we 
found that the surface layer of a cluster is in a melt- 
ing state even at the temperature much lower than the 
melting point of the cluster. The surface melting state is 
almost equivalent to a condition where atoms keep rear- 
ranging along the surface of the cluster. Such an active 
surface motion is converted into the rapid diffusion of 
solute atoms in the direction perpendicular to the sur- 
face and results in a rapid SA. As far as such a diffu- 
sive motion assisted by surface melting is concerned, the 
rapid diffusion of the solute atoms into the solid clus- 
ter is quite different from the diffusion process in a bulk 
solid. It should be emphasized that the surface melting 
is the very origin to activate the radial diffusion pro- 
cess. As discussed in Sec. VI, the radial diffusion process 
is accerelated by a successive gliding motion of surface 
atoms, even if core part of a cluster is solid-like. The ac- 
tive motion of surface atoms is gradually converted into 
the rapid radial diffusion by the frequent onset of a glid- 
ing motion of the surface atoms. A gliding motion as 
a collective atomic motion will be elaborated in detail 
elsewhere p4[ . 

Before closing our conclusion, it is worth recalling three 
important factors, say substrate effect, dimensionality ef- 
fect, and manybody effect, which are not taken into ac- 
count in the present study. In fact we neglect the role of 
substrate which support a cluster and absorbs the heat 
accumulated in the alloying process. By choosing isoe- 
negetic MD simulation we assumed that the coupling 
strength between substrate and a cluster is very week 
and heat transfer from cluster to substrate is consider- 
ably slow. Although we roughly evaluate how fast is the 
heat transport in the appendix, there still exists a pos- 
sibility that we underestimate the effect of a supporting 
substrate. If we emphasize a role of a substrate as a 
heat reservoir for clusters, an isothermal dynamics such 
as a Langevin simulation should be employed to trace 
the time evolution. According to the Langevin dynam- 
ics, kinetic temperature of a cluster does not increase 
as SA proceeds, because the released heat of solution is 
quickly absorbed by the substrate. Thus, we are able to 
remove considerable temperature rise caused by negative 
heat of solution and explore the effect of temperature as 
purely as possible, distinguishing from the effect of nega- 
tive heat of solution [^5| . As far as conclusions we present 
here concerned, the gross feature of SA in an isoenergetic 
condition is not much different from that in an isothermal 
one. 

In the MD studies of bulk metals, manybody poten- 
tial models are usualy employed to mimic interaction be- 
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tween metal atoms. We examined an isoenergetic MD 
of SA with a manybody potential based upon the EAM, 
but we did not observe very significant differences from 
the present simulation p9|] . From these facts we expect 
that, except for some detailed apects, the peculiar fea- 
tures due to manybody potential do not essentially alter 
our results related to dynamics of SA process. In this 
connection, it is worth noting that the presence of spon- 
taneous mixing behavior has been also reported for alkali 
halide microclusters(KBr-KCl system) J2(|. The atomic 
interaction in alkali halide compound is well-described by 
Born-Meyer type pairwise potential, which is completely 
different from manybody interaction of metal atoms [^7| . 
For these reasons it is plausible to say that manybody 
effect, which is peculiar to metal, is not essential for the 
onset of SA. We employed a 2D model which is somehow 
special in a sense that it exhibits anomalous fluctuat- 
ing properties near the melting point, which is similar 
to Kosterlitz-Thouless type transition p8|| . The direct 
outcome due to the confinement in 2D, not 3D, is also 
reported in comparison to the 3D model with EAM [^9| . 
Nevertheless, our preliminary results reveal that there 
are no significant differences between 2D Morse model 
and 3D EAM model, as far as the materials we examined 
here are concerned. 
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APPENDIX A: HOW FAST DOES KINETIC 
ENERGY OF ATOMS TRANSFER FROM A 
CLUSTER TO THE SUBSTRATE? 



In the present paper we assume that the heat given 
by formation of bonding between a deposited impurity 
atom and cluster was released to the substrate slowly. In 
the appendix we evaluate how fast kinetic energy trans- 
fers from atoms in cluster to the substrate in terms of 
a simple one dimensional model. As depicted in Fig.A.l 
the cluster atoms are bounded to the substrate atom via 
harmonic potential, where mass of cluster atoms and sub- 
strate atoms are M and to, respectively. 




substrate cluster 



Fig.A.l A schematic picture of a substrate atom(whitc 
circle) interacting with a cluster atom (shaded circle) via 
harmonic potential. 

Spring constant of harmonic potential between a clus- 
ter and a substrate atom is denoted by K. Atomic fre- 
quencies of a cluster atom and a substrate atom are loq 
and wi, respectively. Langevin equation for two atoms 
locating interface between a substrate and a cluster is 
given by 

xi = -u\x\ + — (x 2 - Xl) (A.l) 
K 

x 2 = -u> Q x 2 - 70*2 + f(t) H (xi-x 2 ) {A.2) 

TO 

where x\ and x 2 are positions of a substrate atoms and 
a cluster atom, 70 is a friction constant, f(t) is an ex- 
ternal random noise, to and M are mass of a substrate 
atom and a cluster atoms, ujq and lo\ are vibrational fre- 
quencies of atoms inside the substrate and the cluster. 
Laplace transform of these two equations is expressed as 

z 2 X 1 - Z x{ 0) - X[ 0) = -a£STi + £(X 2 -X 1 ) (A3) 
and 

z 2 X 2 - - if } = 
-ulX 2 - l0 zX 2 + 7o Af } + F(z) + —(X 1 - X 2 ) (A4), 

TO 

where Xi(z) = C[xi(t)] — J °° x\ (t) exp(— zt)dt, X 2 (z) = 

- / °°a; 2 (t)exp(-zt)di, F(z) = C[f(t)}, X[ a) = 

xi(0) and X^ = x 2 (0). They lead us to the following 
expression for X 2 . 

X 2 = - s+H^ 

m(z 2 +^ + ^+ la z) 

where 

Y (0) 

F(z)=F(z) + (z+ l0 )xW+ 2 . (A6) 

z +UJ + 7k+ ^ z 
Our goal is to give a closed form to evaluate an effective 
friction constant for the variable X\ . The substitution of 
(A. 5) to (A. 3) and its inverse Laplace transform yields, 

xi(t) = -ufx! - + ^ J 6{t-t')xi{t')dt' + g{t), 

(A.7) 
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where 



K 



7o a 



0(t) = ^- exp(— £t)ain(wi) 



w o + - - 4 

m 4 



if ^ 



(A8) 



(A9) 



(A10) 



By introducing a new variable $(f) = J" °° 6(t)cIt, the 
third term of (A7) is rewritten as, 

f 0(t-t') Xl {t')dt' = mt-t') Xl {t%- [ ${t-t') Xl (t')dt'. 

Jo Jo 
Then, the resulting expression for x\ is given by 



(1-$ )Z1--^: f $(T-t') Xl dT + g(t), 

(All) 



Xl = — LO x X\ - 



where 



AT 
M 



and $o = ®{t = 0). When comparing time scale of Q(t) 
to that of xi, it is easy to note that xi oscillates with fre- 
quency u>i and that temporal behavior of $ is dominated 
by a frequency to. These values are determined by the 
frequencies of the substrate atoms and the cluster atoms, 
respectively. The value of u>\ is about 0.2 x 10 14 [sec -1 ] 
for Au, while Q is roughly estimated as 3 x 10 14 [sec -1 ] for 
carbon graphite. We may say characteristic time scale for 
$ and ii is well separated. As a result, the third term 
of rhs of eq.(A.ll) is simplified by extracting x\ out of 
integral. 

On the other hand, the explicit form for <E> is 

K e 2 1 ry n 
$(*) = ^(uicosuit + — smut), (A13) 

muj K + £) 2 

and the so-called frequency shift, say $o, is 



$n = 



K 



m (7o 2 + £) 



(A14) 



In addition, if we assume to hold symmetric relation for 
$, namely $(i) = <j>(— t), and to extend upper limit of 
integral region respect to r from t to oo, then J Q <I>(i — 
r)xidr = xi J °° $(r)dr. If we put /? = J °° $(r)dr, then 
one can get the following expression from eqs.(A9) and 
(A13), 



/3 = 



K 



lo 



(A15) 



By taking into account the relation 70 ~ u>o and Wq > ^, 
it is possible to give the following relation, (3 ~ fj^s- 
Consequently, eq.(A.ll) can be rewritten as 



K 



K 



X\ = —lu iXi — — x\ — -77/3x1 + g(t) 



M 



M' 



(A16) 



Then, we finally obtain a simple expression of the effec- 
tive damping factor 7, 



K n K K 1 



AT 



M m 1 



(A17) 



It is safe to say that the damping factor 7 is small enough, 



since the ratio -2- is estimated as, 



— — 

J_ = ^M ^ Vm ) 2 ( tJl ) 
CJl CJi CJq ^0 



(A18) 



Due to the frequency mismatch at the interface between 
the cluster atoms and the substrate atoms, one can easily 



show the following relations, (-^f-) 2 ~ j^, ( 
and (— ) ~ jfi. Finally, we obtain - 1 



~ \2 

10 ' V UJQ I 

10 -3 



10' 



.,, . - , ,.. , indicating 

that an energy transfer from a cluster atom to a substrate 
is sufficiently slow comparing to time scale of the atomic 
frequency of cluster atoms. 
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